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Abstract. We investigate the dynamics of two-component Bose-Einstein condensates 
(BECs), composed of atoms in two distinct hyperfine states, which are linearly 
coupled by two-photon Raman transitions. The condensate is loaded into a double- 
well potential (DWP). A variety of dynamical behaviors, ranging from regular 
Josephson oscillations, to mixed Rabi-Josephson oscillations and to regimes featuring 
an increasing complexity, are described in terms of a reduced Hamiltonian system with 
four degrees of freedoms, which are the numbers of atoms in each component in the 
left and right potential wells, whose canonically conjugate variables are phases of the 
corresponding wave functions. Using the system with the four degrees of freedom, 
we study the dynamics of fractional imbalances of the two bosonic components, and 
compare the results to direct simulations of the Gross-Pitaevskii equations (GPEs) 
describing the bosonic mixture. We perform this analysis when the fractional imbalance 
oscillates around a zero-time averaged value and in the self-trapping regime as well. 
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1. Introduction 

The study of Josephson oscillations and self-trapping both with a single bosonic 
component [1] and in bosonic binary mixtures [21 |3l HI O El [71 [10] trapped in double- 
well potentials (DWPs) has attracted much interest in the context of the current work 
on ultracold quantum gases. A specific ramification of this topic corresponds to the 
situation in which the two components of the mixture are different hyperfine states of 
the same bosonic atom [21 |3l [6], which may be linearly coupled by an external resonant 
field [11]. This setting suggests a possibility to study the interplay between Josephson 
and Rabi oscillations, the latter being induced by the linear interconversion between the 
components [121 [l3]. This subject was considered in several earlier works [HI [151 [T7] . 
In particular, in Ref. [H] the authors analyzed a crossover between the Josephson and 
Rabi dynamics, using a nonstationary model (with a linearly growing magnetic field), 
which, in the Josephson limit, was reduced to a system of two degrees of freedom. In 
Ref. [I5], a two-degrees-of-freedom model was used too, with the objective of studying 
the quasiparticles' spectrum in the symmetry-broken ground state. As suggested in Ref. 
[16] . there exists a possibility to distinguish between the Rabi and Josephson regimes 
by considering a beam-splitter model based on a nonstationary DWP. An experimental 
implementation of internal bosonic Josephson junctions with Rubidium spinor Bose- 
Einstein condensate has been recently considered by Zibold et al. [18] in connection 
with bifurcations occurring at the transition from Rabi to Josephson dynamics. 

The natural minimum basis for the analysis of the Rabi- Josephson oscillations in 
two-component systems in the DWP should include four degrees of freedom [10]. The 
objective of this work is to analyze physically relevant dynamical regimes within the 
framework of the minimal system. The predictions will be verified via the comparison 
to direct simulations of the underlying Gross-Pitaevskii equations (GPEs). 

The paper is organized as follows. The model is described in Section II, and the 
finite-mode approximation is derived in Section III, where we also verify its accuracy 
by comparison to direct simulations of the underlying GPEs. In Section IV, we report 
the main results of the work, which demonstrate the interplay between the Josephson 
and Rabi oscillations, concluding that the oscillations become characterized by an high 
degree of complexity with the increase of the strength of Rabi coupling. In order to 
enlarge and complete the analysis presented in Ref. [9], the study of the self-trapping 
regime is faced as well. In fact, in Section V - within the self-trapping regime - we 
discuss the accuracy of our model following the same path as in Sec. IV; we comment 
about the influence of the linear coupling constant on the onset of the self-trapping 
dynamics. The paper is concluded by Section VI. 

2. The model 

We consider a binary Bose- Einstein condensate of two different species (with index 1 and 
2) of repulsively interacting bosons. The condensate is trapped in a DWP, which can be 



Rabi- Josephs on oscillations and self-trapped dynamics in atomic junctions with two bosonic species?) 

produced, for example, by a far off-resonance laser barrier separating each component 
into two regions, L (left) and R (right). These components may be , for example, two 
distinct hyperfine states, \F = 2, nip = 1) and \F = l,mF = —1), of ^'^Rb [21 [3]. 
Note that since lAmi^l = 2 one needs two photons to couple the two levels. A weak 
external magnetic field gives rise to a small difference, hwo, between the energy levels 
of these states. Two-photon Raman transitions between the levels, characterized by 
Rabi frequency Q, can be induced by a laser beam of frequency Ud, with detuning 
6 = Ud — ujq. Using the rotating- wave approximation (e.g., neglecting high-frequency 
terms in the atom- field interaction), in the mean- field approximation macroscopic wave 
functions \E'n(r, t) (n = 1, 2) of the two components of the condensates obey the system 
of coupled Gross-Pitaevskii equations (GPEs) fi2\ [T9]: 



where Vtrap(r) is the trapping potential and \i/„(r,t) is subject to the normalization 
condition. 



with Nn{t) the number of bosons of the n-th species. Similarly, m, a„, and (?„ = 
4:7ih^an/mn denote, respectively, the atomic mass, s-wave scattering length, and intra- 
species nonlinearity coefficient of the n-th species (the atomic mass is common for 
both species). The constant accounting for the linear interconversion between the 
bosonic components is expressed in terms of the respective Rabi frequency, Q. Finally, 
gi2 = 27ih'^ai2/m is the coefficient accounting for the nonlinear interaction between the 
species, being the respective s-wave scattering length. Notice that the total number 
of particles Ni(t) + N2{t) is a conserved quantity. In the following, we consider both Qn 
and gi2 as free parameters, due to the possibility to change the scattering lengths by 
means of the Feshbach-resonance technique, see, e.g., Ref. [20] and references therein. 

Equations ([H) can also be derived in a different physical setting, by assuming that 
the two hyperfine states may be coupled by an external ac magnetic field B cos{ujt) of 
frequency u = uq — 6. In this case, the linear coupling term in Eq. ([T]) corresponds to 
the Rabi frequency Q = fi-B/h, where n is the dipole matrix element for the transition 
between the two hyperfine states [TT] . 

The trapping potential for both components is taken to be the superposition of a 
strong harmonic confinement in the transverse {x,y) plane and of a DWP in the axial 
(z) direction, i.e.. 




+ ^12|^3-nP]^n + Y% 



(1) 




(2) 




(3) 



We proceed by writing the Lagrangian associated to the GPEs ([T]), 
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- (^1^2 + ^1^2) - ^12|^in^2 



(4) 



where \l/„ stands for the complex conjugate of To derive, at first, the ID 

approximation, we adopt the usual ansatz 



'^n{x,y,z,t) 



exp 
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y 



2a' 



fn{z,t) 



(5) 



where a±^n = ^/f^/ {jnun) are the respective transverse-confinement radii, with the ID 
wave functions fn{z,t) obeying normalization conditions J dz\fn{z,t)\'' = N„it). Note 
that the factorized ansatz ([5]) is valid under the strong transverse confinement, viz., when 
fl'n|/nP/4vra^ „ <C hun [2I]- By inserting the ansatz ([5]) into the Lagrangian (jlj) and 
performing the integration in the transverse plane, we derive the effective Lagrangian 
for the ID wave functions: 



L,n=/ci.([5:/„(^n| + A.|,)/, 

n=l,2 

{-iYh5 



en + Vbwp(2) 



\fn?-^-^\fn\' 



h/2) (/l/2 + /l/2)-^12|/in/2| 

where the following constants are introduced: 

en = (1/2) [nV {mal^j_) + nnuj^a^^ , 

9n = 9nl (27ra^ „) , 

VL = 2aj_^ia^^2('^i,i + 0^,2)"^^ > 

and 

912 = (^12/vr) {al l + (^±,2y^ ■ 
By varying Lid with respect to /„, we derive the effective ID GPEs, 



(6) 
(7) 
(8) 

(9) 



■' + [en + Vuwp{z) + 



2m„ dz'^ 



+ 9n\fn\' + ~gi2\f3-m\Vn+ h 



(10) 



3. The finite-mode system 



To approximate the dynamics by a finite-mode approximation, we make use of the 
two-mode decomposition for each wave function, as originally introduced in Ref. [22]: 

fn{z,t) = ij^{t)<P^{z)+^P^it)<P^{z) . (11) 
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The orthonormal real functions (p^i^) constructed according to the same path as in 
Ref. [7] and as commented below. These functions are localized in the left and in the 
right wells, respectively {a = L,R) [7j, and 

= v/iV^e^'^-W, (12) 

with the total number of particles in the n-th species being N^{t) + N^(t) = |^ + 

1-0^(^)1^ = Nn(t). As described in the Appendix, we derive explicit evolution equations 
for the temporal evolution of the fractional imbalances, Zn = {N^ — N^)/Nn, and intra- 
species relative phases, 9n = 0^ — 6',^: 
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_„ cos6'3_„ + Kc,i2)Ns_n^/l " sm 9r. 
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Zn COS 26n 



h 

Ul2 - V,2 



h 



COS d 
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3— n 



K. 



c,12 
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AT \ \/ 1 , ■COS7L-W- — COS7ijl . (13) 

2h\l Nn \\l 1 + Zn y 1~ Zn / 

Here the signs plus and minus pertain to = 1 and 2, respectively. The total numbers of 
particles Nn of each component and the respective phases, ■ja = 6*1 —6*2 (recall a = L,R) 
evolve according to: 

Nn = ±(-^ VA^iA^2)(V(1 + -i)(l + Z2) sin7i + ^(1 - 2;i)(l - ^2) sin7R), 



1 
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where AE = E2 — Ei = h6. The latter relation follows from the normalization of 
functions (pa and the first equation from Eq. f l20|) . Obviously, in the absence of Rabi 
coupling, i.e. ^2 = 0, the total number of particles is conserved in each component. 
In that case, Eqs. (fT3l) reduce to the equations of motion for coupled pendula derived 
in Ref. [7]. We observe that due to the nonlinearity associated to the intra- and 
inter-species interactions, the system is nonintegrable also when 1] = 0. On the other 
hand, in the absence of the Josephson coupling, i.e., when Kn = 0, the dynamics will 
be characterized by independent Rabi oscillations in each well; when the Josephson 
coupling is finite, i.e., 7^ 0, and much smaller than Cl, the aforementioned single- well 
independent Rabi oscillations will be weakly coupled by the Josephson tunneling. Note 
that both when Kn = and Kn 7^ 0, the Rabi interconversions will be deformed by 
nonlinear effects due to the intra- and the inter-species interactions. 
In our calculations the axial double-well potential is given by 



Vdwp{z) — —Uq 



sech 



2fZ-\-ZQ, 



+ sech 



zo- 



A 



with 



Uo = hiu±[l + sech.\'^)y 



A 



(15) 



(16) 



that is the combination of two Poschl- Teller (PT) potentials separated by a potential 
barrier the height of which may be changed by changing A, and centered around the 
points —zq and Zq. The the wave functions of the eigenvalues problem in the presence of 
the only potential Va{z) {a = L, R) are exactly known. In particular, the wave function 



(14) 
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of the ground state is 

c^r^{z) = B[l-Tanh\'-^t-'' 

If a = L (R), the function (j)'^'^'^{z) is centered around —zq (+20)5 and it is the ground 
state wave function of the PT potential centered around the point —Zq {zq). In Eq. 
(fTTl) B, equal for both sides, ensures the normalization of the wave function. The above 
functions (p^i^) and (p^i.^)- Then, proceeding from the functions f|T7|) . (pL^z) and (pR^z) 
can be determined following the same perturbative approach as in [7j, where it is shown 
that, under certain hypothesis, the aforementioned functions can be written in terms of 
(pi^^^{z) and0^'^^(^) given by ([I7D. We get [7] 



1 



1 > L,PT/\ I / 1 1 \iR,PT/ 



(18) 

where. = /+~rfz0^'^^(z)0r^(z). 

The relations between the macroscopic parameters involved at right hand sides of 
Eqs. ( fT3l) -( fT4l) and the microscopic parameters of the problem are reported in the Eq. 
(^U^ of the Appendix. We observe, moreover, that to make the system fully symmetric 
we also assume that Ki = K2 = K , Ui = U2 = U , K^^i = _K'c,2 = Kc, and Vi = V2 = V. 

Obviously, in the absence of Rabi coupling, i.e., ^2 = 0, the total number of particles 
is conserved in each component. In that case, Eqs. f|T3l) reduce to the equations of 
motion for coupled pendula derived in Ref. [7]. Note, moreover, that in our calculations 
we have assumed that the both the components feel the same harmonic potential so 
that ui = U2 = uj±. 

To verify the reliability of the finite-mode approximation which leads to the ODEs 
f lT3|) and f|T^ )-f lT^ . we compare the evolution of fractional imbalances Zn, as predicted 
by this system, with results of direct numerical simulations of the ID GPEs ( ITOj) both 
when the fractional imbalances Zn{t) oscillate around a zero-time averaged value and 
when the time-averaged value of {zn(t)) 7^ 0, that is the self-trapping regime. For the 
oscillations characterized by {zn{t)) = 0, the results of the comparison are presented 
in Fig. [TJ This figure shows a good agreement, especially when the Rabi coupling, Cl, 
is small enough. At larger values of ^2, the finite-mode approximation demonstrates a 
deviation, which accumulates at sufficiently long times. 

In Fig. 121 we report the above comparison between GPEs and ODEs when the 
fractional imbalances are both self-trapped. Also in this case the distance between the 
predictions from the two approaches increases for long times when Cl is big enough. It is 
worth to observe that the authors of Ref. [9] have performed the comparison between the 
predictions of GPEs and those ones ODEs deriving from the two-mode approximation 
only when the fractional imbalances oscillate around zero. Julia-Diaz and co-workers. 
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Figure 1. Fractional imbalances zi and Z2 of the two components vs time t in the 
Josephson regime. Solid lines: ID GPE, Eqs. (fTOl) . Dotted lines: finite-mode equations 
(fT3ll - (fT4|) . Upper panels: = (no Rabi coupling); central panels: D, = K; third 
row: il = 3K; lower panels: fl = 2QK. The parameters of the double-well potential 
^ arc A = 1 and zq = 3. Ki = A's = K = 4.955 x lO'^, C/i = C/2 = C/ = O.IK, 
and C/12 = -O.OIC/. Other parameters: = -3.684 x lO"*^, V = 2.268 x 10"^ 
Kc,i2 = -O.Olifc, 1^12 = -0.005T/; E2-Ei = -0.39, as in Ref. [12]. Initial conditions: 
iVi'(0) = 200, iV2(0) = 100 zi(0) = 0.5 = -22(0), 6i„(0) = 0, 7„(0) = 0. Time in units 



of ( 



lengths in units of a±, and energies in units of 



moreover, have integrated the aforementioned ODEs under the assumptions of small 
imbalances, and small intra- and inter-species phase differences [9]. 

4. Josephson and Rabi oscillations mixing and the self-trapped dynamics 

One of the central theme of this work is the interplay of the Josephson oscillations and 
Rabi oscillations. We start with the situation in which only the Josephson coupling is 
present, i.e., = 0, see the upper panels of Fig. [Hand of Figs. [3], HI In the other panels 
of these figures we report the behavior of the system in the presence of finite values of VL. 
One can see that, the greater is the Rabi coupling between the two species (hyperfine 
atomic states), the greater is the dynamical complexity exhibited by the system, as 
observed, especially, in the lower panels of Figs. [Tj [3] and HI Moreover, as shown in the 
lower panels of Fig. [1], when VL = 20K, the behaviors of zi{t) and Z2{t) are strongly 
asymmetric with respect to each other; this asymmetry is absent when E2 — Ei = 0. 
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Figure 2. Fractional imbalances zi and Z2 of the two components vs time t in the self- 
trapping regime. Solid lines: ID GPE, Eqs. (|10l) . Dotted lines: finite- mode equations 
(fT3|) - ITil) . Upper panels: 17 = (no Rabi coupling); central panels: Ct = 5K; lower 
panels: il ^ lOK. We choose Ki = K2 = K = 4.955 x IQ-^, Ui = U2 = U ^ O.IK, 
and U12 ~ — 2C/, A'c.12 — —"^Kc, V12 — —V. Other parameters, initial conditions and 
units as in Fig. [T] 

In Fig. |3] we plot the dynamics of the finite-dimensional system in the planes of 
fractional imbalances Znify and phases 6'„(t), using the trapping and input parameters of 
Fig. [TJ The figure shows the motion is fully periodic for = 0, than it becomes quasi- 
periodic and finally aperiodic by increasing the value of fi. To quantify the increasing 
of complexity with the growth of f2, in Fig. |4] we show the power spectrum of the 
oscillations, represented by the absolute value, |2;„(a;)|, of the Fourier transform of z^it^, 
denoting by Gj\ the frequency associated to the maximum of |2;i(a;)|, and by 0)2 the 
one associated to the maximum of 1 2:2(0;)!. The frequencies GJn are the fundamental 
frequencies of z^it), i.e. the frequencies of the carrier waves of Znif)- Let us focus on 
Table 1 obtained with A'^i(O) = 200 > N^^^ = 100. From the data there reported, we 
see that - whatever is the value of the Rabi coupling Q - ui > UJ2. Then the multi- 
particle tunneling period associated to cui, Ti, given by 2ti/(jJi, is smaller than the one 
associated to ^2-, T2, given by 271/002- This is due to the fact that - within a mechanical 
analogy - the ODEs equations ( fT3i) - (|T4l) describes two coupled pendula. At least in the 
presence of sufficiently weak inter-species interactions, the mass of each pendulum is 
related to the inverse of the particles number Nn of the nth species, as discussed for 
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Figure 3. Josephson regime. The dynamics of the finite-dimensional system in 
the planes of fractional imbalances z„(t) and phases 0„(i). The parameters, initial 
conditions, and units as in Fig. [TJ 



a single component in Ref. [T]. From the Table 1, it is possible infer as well that the 
greater is Q the greater is the ratio T2/T1. On the other hand, if we keep fixed the Rabi 
coupling and increase A'^i(O) with respect to A''2(0), the period Ti will be smaller and 
smaller if compared with T2. Moreover, from Table 1 again, we see that the relative 
changes r„ of Un with respect to their values at = increase by increasing Q. 

From Fig. |4] it is possible to gain physical insight in the dynamics of the system, 
especially for large values of the Rabi coupling Q. We can see that when Q is sufficiently 
large, |2;n(i^)| exhibits a multi-peak structure related to the appearance of frequencies 
different from the fundamental one. This reflects in an increasing number of harmonics 
involved in Zn{t) (see, in particular, the lower panels of Fig. [1] where Q = 20K) and, 
accordingly, in an increasing degree of complexity. By analyzing Fig. H] again, one can 
conclude that the power spectrum approaches that of random oscillations as Q increases. 
This, immediately, reflects on a quite complicated dynamics when we focus on the plane 
Zn(t) — 9n(t) (see, in particular, the lower panels in Fig. [3] where Q = 20K). Finally, 
note that the complexity pertaining to high values of Q increases if the dynamics is 
observed on sufficiently long time scales as shown in the lower panels of Fig. [3] and Fig. 
il 
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Figure 4. Josephson regime. The absolute value of |z„(w)| of the Fourier transform of 
the fractional imbalances vs. frequency uj. Note that the vertical axis is in logarithmic 
scale. The parameters, initial conditions, and units as in Fig. [1] 
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Table 1. Spectra of Fig. HI Frequencies uin of the maxima of for different values of O 

and their relative changes r„, with respect to the absence of the Rabi coupling. 

For the investigation of the self-trapped dynamics let us consider again Fig. |2] and 
also Figs. |5] and [6l In Fig. [5] we report the dynamics in the planes {zn,0), while in 
Fig. E]) we show the absolute value |2;n(a;)| of the Fourier transform of Zn{t). From these 
figures, we see that in correspondence to high value of the linear coupling constant, an 
increasing complexity is observed within the dynamics supported by the junction. 

It is important to stress that it exists a crossover value of the intra-species 
interaction, say Ust, (note that we are setting Ui = U2 = U), such that the two 
components are both self-trapped when U > Ust- By integrating the ODEs (fT3l) - 
f|T^ . we study f/s^ as a function of the absolute value of the inter-species interaction 
amplitude U12 for different values of Q. The results of this study are presented in Fig. 
[71 /^From the plots there reported, we can see that, once fixed U12, the greater is Q, the 
greater is the value of the intra-species interaction necessary to achieve the self-trapping. 
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Figure 5. Self-trapping regime. The dynamies of the finite-dimensional system in 
the planes of fractional imbalances and phases 9n{t). The parameters, initial 

conditions, and units as in Fig. [2l 



5. Conclusion 



In this work, we liave introduced the model which allows one to study the interplay of 
the Josephson and Rabi oscillations in a binary Bose-Einstein condensate trapped in the 
double-well potential structure. The Rabi coupling is provided by the interconversion 
between the two species of the condensate, which represent distinct hyperfine states of 
the same atom. 

To capture core features of the dynamics, we have derived a finite-mode 
approximation with four degrees of freedom, that represent the populations of the two 
species in the two symmetric potential wells. Comparison to full simulations of the 
underlying Gross-Pitaevskii system demonstrates that the truncated system provided 
for a reasonable accuracy. Systematic simulations of the system reveal the transition 
from regular Josephson oscillations to complex dynamics with the increase of the 
Rabi interconversion rate. Within the framework of this analysis, we have discussed 
the possibility to infer, at least at a qualitative level, the behavior of the fractional 
imbalances and the dynamics in the plane Zn{t) — 6n(t) proceeding from the Fourier 
analysis of the fractional imbalances, especially for large Rabi couplings. 

We have analyzed the dynamics of the atomic Josephson junction when both the 
components are self-trapped as well. We have shown that also in this case the truncated 
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Figure 6. Self-trapping regime. The absolute value of of the Fourier transform 

of the fractional imbalances vs. frequency to in the self-trapping regime. Note that the 
vertical axis is in logarithmic scale. The parameters, initial conditions, and units as in 
Fig. El 



system gives rise to reliable predictions since the good agreement with the predictions 
of the associated Gross-Pitaevskii system. We have studied the influence of the linear 
coupling on the self-trapping onset. 



Appendix 

In this appendix, we discuss the path followed in deriving the evolution equations for 
the fractional imbalances Zn = {N^ — N^)/Nn and the intra-species relative phases 
On = On - i-^- Eq. and the Eq. ([H]). We start by deriving the effective 

Lagrangian L^s in terms of variables A^^ and : 



-eff - Yl 

n=l,2 



[ - hOX - MX - EX - E^N^ 



+ 2K^^/NXcos{0^-0^^) 

-2K,,nNn^/N^cosiO^-0^) 
-VnK';K^i2 + cos 2(0^^ -0^^ 



- 2 i^,,i2 {Nf^ + N^) \ N^N^' cos(^f - 0^ 



Rabi-Josephson oscillations and self-trapped dynamics in atomic junctions with two bosonic specieslA 

0.002 



^0.0014 




0.0008 



0.0001 0.0002 

1U12I 



0.0003 



Figure 7. Self-trapping crossover value Ust as a function of the inter-species 
interaction \Ui2\ for different values of (l. Dotted line: 17 = 0; dashed line: f2 = 5K; 
dot-dashed line: (l = lOK. The double-well potential, initial conditions and units as 
in Fig. H 



2{R^,JN[Ni^ cos( 



- 2(/2f2 JiVfiV2« cos( 



qR 



12 Y ^'1 ^'2 ^^^K"2 

where the following constants are introduced: 



\2ma'^, 



dz 



+ V^DWP(^) + 



U: = ~gn dz 



(19) 



Rabi-Josephson oscillations and self-trapped dynamics in atomic junctions with two bosonic specieslb 



K, 



dz 



lL jaR 



2m dz dz 

dz{rn{z)f€{^) 



c,n 9n 

K = 2^.y dzmz)f{<pp^{z)f 

u^2 = ~9i2 1 dz{<P'^y{<f>^)' 

Kc,l2 = ~gi2 j dz{(P'^{z)f4{z) 

v,2 = ~gi2 1 dz{r,{z)f{4{^)f 



^12 

qLR 
^12 



{n/2) / dz<t>1{zW,{z) 



RL 
12 



(fi/2) / dz<Piiz)^!^{z) 



(1/2) I dz(f)f{z)(j)^{z) . 



(20) 



To analyze the finite-mode dynamics induced by Lagrangian (1191) . we define the 
canonical momenta conjugate to generalized coordinates h9'^: 
1 aLeff 



Pes 



(21) 



Accordingly, the Hamiltonian of the system is written in terms of the canonical 
coordinates and momenta as follows: 



n=l,2 

- 5^ 2JCVp^cos(^^-0 



n=l,2 



+ E if + 

n=l,2 

+ K(2 + cos2(^^^-0)PW] 
+ Uf-^PefPe^ + U^2PeRPeR 



2K. 



c,12 



{Pe^ + Pe^)^P0^PeR cos{9^ - 9^ ) 

+ {Pe^ + Pe^)^Pe{Pe^ cos(ef - 9^) 
+ yi2{PeiPeR + PefPe^) 

+ 4 Vi2{^/p^^/n^ cos(^f - 9{) cos(^f - 9^)) 
+ 2iR^,^/p^ cosi9^ - 9^)) 
+ 2(i?f2^/p^cos(^f-^f)) 
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+ 2{T,f^/p^ cos(^2^ - 9f)) . (22) 

The evolution equations for populations A^" and phases 9^ are derived, as the 
canonical equations, from the Hamiltonian: (!22|) 

1 dH ■ 1 dH 



9^ 



h dpeo; 



(23) 



We observe that due to the ortho normality of the decomposition basis, = f^/2 and 
^12' — '^21^ — 0- From the symmetry between the two wells in the DWP structure it 
also follows that = ^ En, = = Un, U^^ = Uf^ = Uu- Using Eq. ([23]), we 
derive explicit evolution equations for the Eq. ( IT3l) and the evolution equations for the 
total numbers of particles of each component Nn and the respective phases, 'ja = 9f — 6*2 
(recall that a = L,R): 



2{Kn - K^^^Nn 



h 



zl sin 9n 



2h 



l-zl)sm29n 



cos6'3„„ + i^c,i2)iV'3-nV 1 " 4 sm 6*^ 



2h 

Ur,. - Vr 



A/iViiV2(V(l + zi){l + Z2)sin7i - a/(1 - zi){l -Z2)sin7R 

2{Kn - Kc,nNn) Z^ COS 9n 



-NnZn + 



h 

Zn COS 29n + 



2h 



h 

Ul2 - Vl2 

h 



Ns-nZS- 



+ 



2 

h 



12 



1 - z|„„ cos 6'3_„ + K, 



c,12 



N- 



Zn COS 9r, 



3— n 



N- 



3— n 



2hV N„ 



1 + Z3-r 
l + Zn 



■ COS 71, - 



2^3- 



cos7ij 



(24) 



iV„ = V^i^2)(\/(1 + zi){l + Z2) sin7L + ^/{l-z{){l-Z2) sin 7^), 

7L = ^(^1(^12 - U){1 + 21) - iV2(t/i2 - f/)(l + Z2) + A^) 



+ — { Ki \ I — cos 9i — Ko \ / cos 90 



n 



iV2(l + Z2)-iVl(l + ^i; 



2hVNj^\ v/(i + ^0(i + ^2) 



COS 7l 



+ ^[iVi^l - cos^i - N2\/l - ^2cos^2] - 



v_ 

2h 
1 



A^i(l - zi){2 + cos 2^1) - iV2(l - Z2){2 + cos 2^2) 



i^e,lA^l(2 + Zi) + i^e,12A^2 + ^12iV2\/l - ^2' COS ^2 



1-^1 
1 + 2:1 



COS 6*1 + 
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Ke,2iV2(2 + ^2)+i^c,12iVl 



1-22 
1 + 22 



COS^o 



+ Vl2Nl^J I - zl cos Oi 

^(iVi(f/i2 - U){1 - zi) - N2{Uu - U){1 - Z2) + Ae) 



1 



A"i 4/ ^— ^ cos 6^1 - K2\ ] cos 6^2 

AlV VI — -21 \ 1 — Z2 



n 



N2{1- Z2)-m{l-Z,] 
v/(l-Zl)(l-22) 



cos7ij 



+ - 2? cos^i - iV2\/l - .zicos^: 



h 

V 

2h 
1 



Ni{l + zi)(2 + cos 2^1) - iV2(l + 22)(2 + cos 2^2 



Ke,lA^l(2 - 2l) + Kc,12iV2 + Vi2N2\ll - zl COS 62 

K,^2N2i2 - Z2) + K^^uNi 



1 + 21 

1-^1 



COS 6'i + 



+ ■^^i2A^i\/l - ^iCOS^i 



i±ii cOSe2 + ^(iVl(l + ;2l) - N2{1 + Z2)), 
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